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ABSTRACT 



Aims. In this paper we study the Spitzer and TIMMI2 infrared spectra of post-AGB disc sources, both in the Galaxy and the LMC. 
Using the observed infrared spectra we determine the mineralogy and dust parameters of the discs, and look for possible differences 
between the Galactic and extragalactic sources. 

Methods. Modelling the full spectral range observed allows us to determine the dust species present in the disc and different physical 
parameters such as grain sizes, dust abundance ratios, and the dust and continuum temperatures. 

Results. We find that all the discs are dominated by emission features of crystalline and amorphous silicate dust. Only a few sample 
sources show features due to CO2 gas or carbonaceous molecules such as PAHs and Ceo fullerenes. Our analysis shows that dust grain 
processing in these discs is strong, resulting in large average grain sizes and a very high crystallinity fraction. However, we do not find 
any correlations between the derived dust parameters and properties of the central source. There also does not seem to be a noticeable 
difference between the mineralogy of the Galactic and LMC sources. Even though the observed spectra are very similar to those of 
protoplanetary discs around young stars, showing similar mineralogy and strong grain processing, we do find evidence for differences 
in the physical and chemical processes of the dust processing. 

Key words, stars: AGB, post-AGB - stars: binaries - stars: circumstellar matter - stars: abundances - Magellanic Clouds 



1. Introduction 

Studies of the chemistry and geometry of circumstellar discs 
have, so far, mainly focussed on the protoplanetary discs around 
young stars (e.g. Meeus et al. 2001; Bouwmanetal. 2008; 



lJuhasz et al.li20 10). However, in recent years it became clear that 
circumstellar discs are present in nearly all stages of stellar evo- 
lution, goingjrom first-ascent giants (Jura 2003; Verhoelst et al. 
I2OO7L iMehs et al.ll2010l) . B[e] supergiants (Kastneret al. 2010), 
asym p totic giant branch (AGB ) stars (e.g. Yamamura et al. 
I2OOOI: IChiu e t al.' '200 61: iDeroo e t al. 2007| ) (proto-)planetary 
nebulae (e.g. Chesneau et al.' '2006, 2007'riLvkou et al 1 1201 ll) 
to white dwarves (e.g. Becklin et al.l 12001 iDong et al. l l2oToh . 
Even though circumstellar discs appear common throughout the 
Hertzsprung-Russell diagram, it is still unclear what links the 
different disc-bearing objects throughout all the late evolution- 
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ary stages. It is likely that there are different formation channels 
depending on the evolutionary status of the central object. 

Whereas for young stars the disc is a by-product of the star 
formation, there is evidence that for the majority of the evolved 
stars the disc is newly formed. The exact formation mechanisms 
are unknown, and will most likely differ for different evolution- 
ary stages. For example, disc formation has been linked to binary 
mergers, wind capture or Roche-lobe overflow (see references 
above). However, in most cases, binarity appears to be the key 
ingredient to the formation of discs in later stages of stellar evo- 
lution. 

In this work we study a particular class of evolved bi- 
nary post-AGB stars surrounded by stable dusty discs. These 
sources were initially selected on the basis of their very strong 
near-infrared excess. Follow-up studies confirmed the bina- 
rity, and showed that the companion star is most likely a 
main-sequence star, with a typical separation of about 1 AU 
dWn Winckel et al ] l2009h . The presence of a disc was already 

S >roposed to explain th e presence of hot dust in the system 
De Ruvte r et al. '2006") and later resolved by interferometric 
observ ations (Buiarrabal et al. 2001, 2007; Deroo et al. 200^ 
Deroo 2007). The discs also explain the observed depletion pro- 
cess in the photospheric abundances of the central post-AGB star 
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(IWaters et al.ll992l:lMaas et al.l2005HGrelen et al.l2009bh . Since 
the dust sublimation radii for these sources are well beyond the 
orbit, all the discs are circumbinary. 

Our previous studies have shown that the discs are ideal envi- 
ronments for strong du st process ing, in the form of grain growth 
and crystallisation (Gi elen et alfc OOS. 2009a b). This dust com- 
position is very similar to what is observed for protoplanetary 
discs around young stars, even though the disc formation mech- 
anisms, and probably also the initial dust species, are very dif- 
ferent. 

In the Galaxy, aro und 80 such systems are now known 
dPe Ruvter et al.l2006h . Recently, large programmes, such as the 
Spitzer SAGE (Surveying the Agents of Galaxy Evolution) pho- 
tometric (Meixner et al. 2006), and follow-up SAGE-Spec spec- 
troscopic (Kempe r et al.ii2010l) . programme indicate that also in 
the Large Magellanic C loud (LMC) pos t -AGB disc sources are 
common: The study of Ivan Aarle et a D (120 llh lists about 650 
probable post-AGB disc candidates in the LMC, and about the 
same number for post-AGB stars surrounded by a cool expand- 
ing dust shell, using SAGE photometric data. 

In this paper we look in more detail to the mineralogy of 
the circumbinary discs, both for sources in the Galaxy and in 
the LMC. For this we use high- and low-resolution Spitzer and 
TIMMI2 infrared spectra. These spectra allow us to study dust 
and gas emission features in the 5-35 fim region. 

The outline of the paper is as follows: In Sections|2]and|3]we 
describe the selected samples Galactic and LMC stars and the 
data reduction process. In Sect. |4] we take a first look at the dif- 
ferent emission features in individual sources, and compare the 
Galactic and LMC sample. The results on the dust parameters 
using a more detailed model to fit the full Spitzer wavelength 
range are described in Sect. |5] Finally, we end with a discussion 
and conclusions in Sects. |6]and Q 



2. Programme Stars 

In this paper we study a total of 57 post-AGB stars with 
evidence for the presence of a stable circumbinary disc, lo- 
cated in the Galaxy and the LMC. The Galactic sample 
consists of 33 stars from the larger sample discussed in 
iDe Ruvter et al.l (l2006l). Of these star s, 21 sources are alr eady 
discussed in iGielen et al] (l2008h and iGielen et al] (l2009ah . To 
complement these 21 sources we obtained Spitzer high- and 
low-resolution spectra of 13 additional suspected post-AGB 
disc sources. The LMC sample consis ts of 24 so u rces, o f 
which 3 are already discussed briefly in IGielen et al.l ( l2009bl) . 
These sources were observed in low-resolution mode , either 
as par t of the larger SAGE-Spec programme of Kempe ret al.l 
(l20Toh . a follow-up to the photometric SAGE legacy pro- 
gramme (iMeixner et al.ll20d6 !). or as part of Spitzer programmes 
3274 (PI: Hans Van Winckel) and 50092 (PI: CUo Gielen) 
(jhttp://irsa.ipac.caltech.edu/data/SPITZER/docs/). 

From .Woods et al., (.201 1.) , we selected the stars which are 
classified as oxygen-rich post-AGB or RV Tauri sources in the 
SAGE-Spec catalogue. We removed the sources for which only 
a small part of the Spitzer wavelength range was observed (LHor 
120-N145 and MACHO 8 1.9728. 14). After this, 16 stars re- 
mained. To increase this LMC sample, we searched the SAGE 
photometric catalogue for the presence of other possible disc 
bearing post-AGB sources. All objects with 24 /vm fluxes be- 
tween 2 mJy and 1 Jy were selected, in order to exclude young 
stellar objects and supergiants. Other selection criteria were cho- 
sen to distinguish between post-AGB stars with an expanding 



shell (F24 > Eg) and binary post-AGB sources with a circumbi- 
nary disc (F24 > 0.5 Eg and J - K < 1). For a detailed description 
of the selection criteria we refer to Ivan A^ le et al. (2011). After 
cross-correlation with optical photometric catalogues and the 
SIMBAD Astronomical Database, 650 sources remained. Of this 
larger sample, the 8 brightest stars were selected and observed 
with the Spitzer infrared spectrograph. For 18 of the 24 LMC 
sources additional ground-based optical spectra were obtained 
at Siding Spring Observatory, the South African Astronomical 
Observatory (SAAO) or with the UVES spectrograph in Paranal. 
This allows us to determine a spectral type and assign an effec- 
tive temperature ( Gielen et a"l]|2009bt Ivan Aarle et"al1l201 1). 

For 15 Galactic sources the binarity has been confirmed 
by radial velocity monitori ng, resulting in orbital p eriods be- 
tween 200 and 1800 days (IVan Winckel et al.1 12009 ). For the 
other Galactic sources, binarity can already be suspected from 
the monitoring programme but not enough data are available to 
derive the exact orbital parameters. Unfortunately, such long- 
term radial velocity monitoring programme for the LMC sources 
is difficult, since it requires several years of observations with a 
high-resolution optical spectrograph on a large telescope, such 
as UVES on the VLT. But, given the strong resemblance of the 
LMC disc candidates to the Galactic disc sources, in chemistry 
of the central star, spectral energy distribution and mineralogy of 
the circumstellar environment (Revniers & Van Winckelll2007l: 
Gielen et al. 2009b), we postulate that these sources will also be 
part of a binary system. 

2. 1 . Spectral energy distribution 

For all Galactic sample stars, spectral energy distributions 
(SEDs) were calculated fr om the photometric da ta and stel- 
lar parameters as given in i De Ruvter et al.l (l2006h . the SAGE 
photometric catalogue and/or the Vizier database. The result- 
ing SEDs can be seen in Figures [T] and |2] From the SED we 
also calculated the luminosity ratio L/r/L,. The total extinction 
E(B - V),oi was determined by dereddening the observed pho- 
tometry and infrared spectra, using the average extinction law 
of Savage & Mathis ( 1 979J extended with the theoretical extinc- 
tion law of lSteenman & Tliel (Il989lll99ll) . Minimising the dif- 
ference between the dereddened observed optical fluxes and the 
appropriate Kurucz model (Kurucz 1979) gives the total colour 
excess E{B - V),ot (Tables[T]|2|. This is done under the assump- 
tion that the extinction is fully due to interstellar extinction, or 
that the circumstellar component follows the same extinction 
law. Since the total extinction probably consists of both an in- 
terstellar and a circumstellar component, the applied deredden- 
ing is thus a maximal correction. The errors on the value for 
E(B - V),ot are calculated using a Monte-Carlo simulation on 
the photometric data. We use an error of 0.05 for the photometric 
measurements in a Gaussian distribution. Since we do not know 
the distances to the Galactic sources, we adopt a likely lumi- 
nosity for evolved low-gravity objects of L, = 5000 + 2000 Lq. 
For the LMC sources we calculate t he luminosity assum ing a 
typical LMC distance of 50000pc (iKemper et all 1201 Ol) . For 
th e LMC sources we use the effective temperatures as given 
in Ivan Aa rle et al. (2011), if available. Since the metallicity and 
\ogg values for these sources are not determined, we used values 
of [Fe/H]= -1.0 and \ogg =1.0 for all stars. These values are 
consistent with those found for the Galactic objects, and have 
only minimal impact on the derived total reddening and infrared 
energy ratio. Of the LMC sample, 6 sources lack optical spectra, 
and thus effective temperatures, and we could not determine the 
total reddening. 



C. Gielen et al.: Silicate features in Galactic and extragalactic post-AGB disks. 



cn 

O 





-7 
-8 
-9 
-10 

-7 
-8 
-9 
-10 

-7 
-8 
-9 
-10 

-7 
-8 
-9 
-10 



-9 
-10 

-7 
-8 
-9 
-10 

-7 
-8 
-9 
-10 



10 0. 



10 0. 



1 10 0. 



10 0. 



EPLyr 


HD131356 

. i \ 

j 


HD7'3985 

4 \ 


HD52961 

g TO 

1 's: 

« \ 


RAS05208 

?^ X 

■ 1 \ ' 




RAS06034 
i" \ o 


RAS06072 

i \ " 




RAS06338 

9 

,v 

'A...^ 


RAS09060 


RAS091 44 

■ J ; 


RAS09538 

i" \ 1 

tl \ ' 

.11 ;. 


|||IRAS1017'1 

O ^ " 


f^. IRAS' 1 000 


IRAS^3258 

o ^ 
O o " 

o 
o 


RAS' 5556 


RASl 6230 


IRAS1 7038 

f \ ^ ' 

r '\ 


RAS' 7233 

o ^ 


IRAS'72^3 

,<^^ 

■ 1 


RASl 7530 

r , N 


RASl 81 23 
;( 


IRAS181 58 
ji' 


IRAS".9125 


RASl 9 157 


RAS20056 

O ""'^-v 

i,« <> 

i' N 






SAO' 73329 


STPjp 


SUCer-i 

g^-. oo 


SXCen 


TWCan 

J 


LYAra 


JYCMa 





-7 
-8 
-9 
-10 



-10 

-7 
-8 



-10 

-7 
-8 
-9 
-10 

-7 

-8 

-9 

-10 

-7 



Wavelength [/xm] 

Fig. 1. The spectral energy distributions of our sample stars. The dereddened fluxes (diamonds), reddened fluxes (gray triangles) 
and Spitzer spectra (solid line) are given together with the scaled photospheric Kurucz model (dashed line). For the sources where 
we lack the stellar parameters to determine the underlying Kurucz model, we only plot the reddened data. 



For some sourc es there is evidence that the d iscs are seen 
close to edge on ( Menzies & Whitelockl 119881: iLlovd EvansI 
11997.) . The visible light of the central source is then seen in 
reflection which makes an accurate determination of the total 
extinction, the luminosity ratio, and the distance very difficult. 
These sources are marked with an asterisk in Tabled] 



3. Observations and data reduction 

3.1. Spitzer 

The spectra were obtained using the SL (A - 5.3 - 14.5 //m), 
LL (^ = 14 - 38 /^m), SH (A = 9.9 - 19.5 //m) and LH 
A = 19.3 - 37 jum) staring modes on th e Spitzer-IRS instrument 
Werner et al.l2004tlHo"uck et al.l2004l) . For the Galactic objects, 
exposure times were chosen to achieve a S/N ratio of 400. For 
the extragalactic objects in our own observing proposal, expo- 
sure times were chosen to give a S/N ratio of 100 for the SL 
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Fig. 2. Same as previous figure. 



mode and 20 for the LL mode. The SAGE-Spec LMC objects 
have a S/N ratio ~ 60 in SL mode and ~ 30 in LL mode. 

The newly obtained spectra from our own Spitzer observa- 
tions were extracted from the SSC raw data pipeline S18.0 ver- 
sion products, using the c2d and feps data reduction packages. 
F or a detailed descrip tion of these reduction packages, we refer 
to lLahuis et al. I (l2006h and lHines et all Jlool . The reduction in- 
cludes background and bad-pixel correction, extraction, defring- 
ing and order matching. Individual orders are corrected for off- 
sets, if necessary, by applying small scaling corrections to match 
the bluer order. For a detailed description of the target selection, 
observi ng strategy and reduc tion of the SAGE-Spec objects, we 
refer to iKemper etal] (|20 1 Oh . 

3.2. TIMMI2 

For some stars we lack the Spitzer IRS-SH observations and 
we obtained additional ground-based N-band infrared spectra 
with the Thermal Infrared Multi Mode Instrument 2 (TIMMI2, 



iReimann etal] 120001: iKaufl et al.ll2003h . mounted on the 3.6m 
telescope at the ESO La Silla Observatory. The low-resolution 
(R ~ 160) N band grism was used in combination with a 
1 .2 arcsec slit; the pixel scale in the spectroscopic mode of 
TIMMI2 is 0.45 arcsec. For the reduction of the spectra we used 
the method described in Ivan Boekel et all (f2005 ). We scaled the 
TIMMI2 spectra to the Spitzer spectra and found a very good 
agreement in spectral shape between the two data sets. 

The resulting spectra can be found in Figure[3]and Figs. IA.6l - 
lAjOl 



4. First inspection of emission features 

Looking at the spectra of the Galactic and LMC sources 
(Figs. IA.6IIA.10b . we find that all sources show clear silicate 
emission. For nearly all sources the prominent broad amorphous 
silicate features at 10 and 20 /im stand out. Furthermore, most 
spectra show additional narrower features at 11.3,6, 19,23,27 
and 33 fim, which are due to crystalline silicate emission. 
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Table 1. The name, equatorial coordinates a and 6 (J2000), effectiv e temperature Teg , surfa ce gravity \ogg and metallicity [Fe/H] 
of our Galactic sample stars. For the model parameters we refer to De Ruvter et al.l (12006 '). Also given is the orbital period (see 
references in De Ruvter et al. 2006; Gielen et al. 2007; Van Winckel et al. 2009). The total reddening E(B - V)tot, the energy ratio 
Liit/Lt and the calculated distance, assuming a luminosity of L» = 5000 + 2000 L©. Stars marked with * are seen in reflection only, 
resulting in unreliable E(B - V),ot values and luminosity ratios, and upper limits for the distances. The last column lists whether the 
spectra are part of the SAGE-Spec catalogue, Spitzer programme 3274 or 50092. 
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-34 44 32 


6500 


1.5 


0.0 




0.78±0.03 


13+9 


10+2.3 


50092 


23 


IRAS 19125+0343 


19 15 00.8 


+03 48 41 


7750 


1.0 


-0.5 


517 


1.08±0.02 


52+3 


1.8+04 


3274 


24 


IRAS 19157-0247 


19 18 22.5 


-02 42 09 


7750 


1.0 


0.0 


120.5 


0.68±0.01 


63+2 


4.2+0.9 


3274 


25 


IRAS 20056+1834* 


20 07 54.8 


+ 18 42 57 


5850 


0.7 


-0.4 




0.51±0.02 


905+42 


10.9+2.3 


3274 


26 


RUCen 


12 09 23.7 


-45 25 35 


6000 


1.5 


-2.0 


1489 


0.55+0.01 


13+1 


2.3+0.5 


3274 


27 


SAO 173329 


07 16 08.3 


-23 27 02 


7000 


1.5 


-0.8 


115.9 


0.39+0.01 


36+1 


6.5+1.3 


3274 


28 


ST Pup 


06 48 56.4 


-37 16 33 


5750 


0.5 


-1.5 


410 


0.00±0.00 


55+1 


5.7+1.2 


3274 


29 


SUGem* 


06 14 00.8 


+27 42 12 


5750 


1.125 


-0.7 




0.58+0.02 


111+7 


4.8+1.0 


3274 


30 


SXCen 


12 21 12.6 


-49 12 41 


6000 


1.0 


-1.0 


600 


0.32+0.02 


34+2 


3.8+0.7 


3274 


31 


TWCam 


04 20 48.1 


+57 26 26 


4800 


0.0 


-0.5 




0.40+0.02 


42+3 


3.2+0.6 


3274 


32 


UYAra* 


17 29 28.9 


-59 54 02 


5500 


0.5 


-1.0 




0.00+0.00 


72+3 


12+2.5 


50092 


33 


UYCMa* 


06 18 16.4 


-17 02 35 


5500 


1.0 


0.0 




0.00+0.00 


89+3 


9.6+2.0 


3274 



IRAS16230 HV12631 




10 15 20 25 30 35 10 20 30 

\ (yitm) A (/um) 



Fig. 3. Best model fits for two of our sample stars, showing the contribution of the different dust species. Top: The observed spectrum 
(black curve) is plotted together with the best model fit (red curve) and the continuum (black solid line). Forsterite is plotted in green, 
enstatite in blue, silica in cyan and amorphous olivine and pyroxene in magenta. Small grains (0. 1 yum) are plotted as dashed lines 
and larger grains (2 and 4/im) as dotted lines. Bottom: The normalised residuals after subtraction of our best model of the observed 
spectra. The models for the other sample stars can be found in Figs. IA.6HA.10l 



6 



C. Gielen et al.: Silicate features in Galactic and extragalactic post-AGB disks. 



Table 2. The name, e quatorial coordinates a and 6 (J2000), and effective temperature Teff our LMC sample stars, taken from 
(Ivan Aarle et aLlUoi ih . The total reddening E(B - V)iot, the energy ratio Lm/Lt, and the luminosity as calculated from our SED 
modelling. The last column lists whether the spectra are part of the SAGE-Spec catalogue, Spitzer programme 3274 or 50092. 



N° 


Name 


a (J2000) 


S (J2000) 




L, 


EiB - VU 


LmlU 


Prog. ID 






(h m s) 


C ' ") 


(K) 


Lo 




(%) 




34 


HV 12631 


05 39 33.1 


-71 21 55 










SAGE-Spec 


35 


HV 2281 


05 03 05.0 


-68 40 25 


5750 


2000 


0.02±0.02 


63+2 


SAGE-Spec 


36 


HV 2444 


05 18 46.0 


-69 03 22 


6750 


4000 


0.26±0.02 


32+2 


SAGE-Spec 


37 


HV2522 


05 26 27.2 


-66 42 59 


6250 


3700 


0.17±0.02 


45±3 


SAGE-Spec 


38 


HV2862 


05 51 21.1 


-69 53 47 


5750 


2700 


0.09±0.02 


54+2 


SAGE-Spec 


39 


HV5829 


05 25 19.3 


-70 54 07 


5500 


1800 


0.00±0.02 


60±2 


SAGE-Spec 


40 


HV915 


05 14 18.0 


-69 12 35 


6250 


4600 


0.21±0.01 


64±3 


SAGE-Spec 


41 


J044458. 18-703522.8 


04 44 58.4 


-70 35 23 


7000 


1400 


0.04±0.02 


40±2 


50092 


42 


J045242.93-704737.4 


04 52 43.2 


-70 47 37 


5500 


2500 


0.30±0.03 


22±3 


50092 


43 


J050143. 18-694048.7 


05 01 43.5 


-69 40 48 


5000 


2700 


0.22±0.02 


19±2 


50092 


44 


J051159.11-692532.8 


05 11 59.4 


-69 25 33 


6250 


8500 


0.02±0.02 


23±2 


50092 


45 


J051333.74-663419.1 


05 13 33.7 


-66 43 19 


6500 


17000 


0.12±0.03 


18±2 


50092 


46 


J052220.87-655551.6 


05 22 21.1 


-65 55 52 


4250 


5000 


0.16±0.03 


2+2 


50092 


47 


J053605.56-695802.9 


05 36 05.9 


-69 58 03 


6750 


8500 


0.38±0.03 


15±3 


50092 


48 


J054312.52-683356.9 


05 43 12.9 


-68 33 57 


6250 


3000 


0.27±0.02 


49±5 


50092 


49 


MACHO 78.6698.38 


05 21 49.1 


-70 04 34 


7000 


3300 


0.46±0.03 


27±2 


SAGE-Spec 


50 


MACHO 82.8405.15 


05 31 50.9 


-69 1 1 46 


6000 


3600 


0.05±0.01 


84±3 


SAGE-Spec 


51 


MSX949 


05 40 14.8 


-69 28 49 










SAGE-Spec 


52 


NGC 1805 SAGE IRS 1 


05 02 24.2 


-66 06 37 










SAGE-Spec 


53 


SAGE050830 


05 08 30.6 


-69 22 37 










SAGE-Spec 


54 


SAGE051453 


05 14 18.2 


-69 17 24 


4250 


2100 


0.10±0.02 


28±2 


SAGE-Spec 


55 


SAGE052707 


05 27 07.2 


-70 20 02 










SAGE-Spec 


56 


SAGE052747 


05 27 47.6 


-71 48 53 










SAGE-Spec 


57 


SAGE054310 


05 43 10.9 


-67 27 28 


4000 


10500 


0.32±0.02 


27±2 


SAGE-Spec 



Even though all the spectra are dominated by oxygen-rich 
dust species, some stars do show evidence for the presence 
of carbonaceous molecules. Clear PAH emission can be seen 
in EPLyr, IRAS 06338 and IRAS 13258, with peaks at 8 and 
1 \2jjm_/Yh&J'A\{ features of EP Lyr were already discussed 
in iGielen et al.l (|2p09a). The peculiar spectrum of IRAS 06338 
not only shows the typical PAH bands, but several smaller fea- 
tures between 6 and 8 /im, most likely resulting from very small 
PAH grains. In this star, the strong narrow peaks between 13 
and 18;um are due to CO2 gas emission, which can also be 
seen in EPLyr and HD 52961. HD 52961 and IRAS 06338 both 
show a strong feature at 18.7yum, w hich can be identified as 
Cgo fullerene emission dCami et alJ201Q') . The detection of these 
carbonaceous molecules in our sample stars will be further dis- 
cussed in an upcoming paper. 

To study the silicate signatures in the infrared spectra, we di- 
vided the full spectrum into 7 different complexes where strong 
silicate emission is seen, more specifically at 10 - 14 - 16 - 
19 - 23 - 27 and 33 ^m. To compare the Galactic stars to the 
LMC sources, we calculated for each group a mean continuum- 
subtracted spectrum in these 7 complexes. The continuum was 
determined by linearly interpolating between the beginning and 
end of the studied regions. These mean spectra are then nor- 
malised to the maximum flux in the wavelength interval. An 
overview of the different mean spectra can be seen in Figure |4] 
together with synthetic spectra of crystalline and amorphous sil- 
icates. Below we discuss the different complexes in more detail, 
and results can be seen in Figure|6]and Figs. lA.ll to IA.5I 

4.1. The Wfim complex (8 - 13 fim ) 

It is clear from Fig. |6] that the mean spectrum in this region is 
very similar for the Galactic and LMC sources. A flat-topped 
feature is seen, where the two peaks come from the emission 



of amorphous and crystalline olivine. The observed 10 /im com- 
plex seems to be broader at the left shoulder, compared to the 
emission feature of amorphous olivine. Additional emission near 
9 fim could point to the presence of amorphous pyroxene or sil- 
ica, which peak at shorter wavelengths. 

Some individual sources do not follow the calculated mean 
complex. IRAS 13258 and EPLyr show no silicate features, but 
exhibit emission due to PAHs. PAH emission probably also con- 
tributes to features seen in IRAS 06338 and HD 52961. Also 
note that in IRAS 06338, the strong feature at 9yL(m seems to be 
shifted bluewards in comparison to the mean. This could point 
to the dominance of silica in this source. IRAS 10174 shows al- 
most no emission of crystalline species, not only at lOjum but 
along its entire wavelength range, and is very similar to the ex- 
tragalactic source J051333. These two sources also do not show 
the broadening at the left shoulder of the complex, and are thus 
expected to be devoid of sihca. 

The lOfim complex is a good tracer of gra in processing, in 
the form of grain size and crystallisation ( van Boekel et al .I2003L 
120051: lJuhasz et al.ll2010l) . Since amorphous and crystalline sil- 
icates peak at two distinct wavelengths, respectively 9.8 and 
11.3 yum, the continuum-subtracted 1 1 .3/9.8 fim flux can be used 
as a measure for the amount for the crystallisation the dust 
has undergone. Furthermore, the peak-to-continuum ratio of the 
lQfj.m complex can be used as a tracer for grain growth, since 
larger grains will result in a less pronounced feature. In Figure|5] 
we plot these two ratios. We do not plot EP Lyr and IRAS 13258, 
since they show strong PAH emission at 1 1.3 fim, contaminating 
the crystalline emission at this wavelength. SAGE 050830 has a 
very high peak-to-continuum ratio of 5.05 ( with a 11.3/9.8/im 
ratio of 0.81), and falls outside our plot range. Most of our 
sources show rather high 11.3/9.8/im ratios, with low peak-to- 
continuum values, showing that the crystallinity fraction is high, 
and the average grain sizes relatively large. 
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Fig. 4. Overview of the different emission complexes. Top: In black we plot the normalised mean continuum-subtracted Galactic 
spectrum, in red the normalised mean LMC spectrum. We did not include the 14 and 16/im complexes, since the noise level for the 
LMC sources made it impossible to determine a mean spectrum. Bottom: The different normalised continuum-subtracted spectra of 
forsterite, enstatite and amorphous olivine are given in respectively blue, green and magenta. 




2.0 2.5 3.0 

Peak to continuum 

Fig. 5. Ratio of the continuum-subtracted flux at 9.8 and 11.3 /im 
versus the peak-to-continuum ratio of the lO/zm silicate feature. 
Galactic sources are plotted in red plus signs, LMC sources in 
blue diamonds. The gray area shows typical values found for 
protoplanetary discs around young stars. The numbers corre- 
spond to numbers given in Tables [1] and |2] 



A very weak correlation (Kendall rank correlation r - 
-0.23), can be seen. This is in contrast to the strong correlation 
seen in discs around young stars between grain g rowth a nd crys- 
tallisa tion processes (Ivan Boekel et al.ll200ll2005; Juhasz et al.l 
The gray area in Figure |5] shows typical values found 
for protoplanetary discs, and it is clear that our sources show 
a much larger spread in values for the continuum-subtracted 
11.3/9.8/im flux ratio. This could mean that in the case of the 
post-AGB discs, the dust might not consist of very small (0. 1 fim 
) amorphous grains, but may already have a higher crystallinity 
or larger grain size. It could also mean that different grain pro- 
cesses are at play, resulting in a slightly different dust grain evo- 
lution. We find no evidence for different behaviour between our 
Galactic and LMC samples. 



4.2. The 14nm and 16 complexes ( 13.5 - 15jum and 
15 - lli^m) 

This region is dominated by two different emission complexes, 
respectively around 14 and 16 fim , as can be seen in Figure [ATI 
Because of the high noise level in the LMC sources for this re- 
gion, we could only calculate a mean spectrum for the Galactic 
stars. 

The I4fj.m complex is sensitive to the emission of enstatite, 
which shows a clear feature around 13.8jum in our observed 
sources. The predicted feature at 14.4 //m, however, is not seen. 
Instead, we do see a clear signature around 14.7yum. The syn- 
thetic spectra of enstatite are known to be sensitive to the re- 
fractory indic es used and th e adopt ed grain size (s ee for exam- 
ple Fig. 10 in Molster et al.' ('2002a') and Fig. 20 in ' Juhasz etall 
0(2010)). Chihara et al. (2002) present an overview of the shift of 
peak position of crystalline pyroxenes with different iron contri- 
butions and we find that the peak positions found in our spectra 
are better modelled with enstatite with a small iron contribution 
of about 10%. In Figure Q we show the continuum-subtracted 
spectrum of ST Pup, which has the strongest enstatite features 
and best S/N ratio in this region of our sample stars, together 
with the laboratory spectra of ortho-ensta tite and clino-enstatit e 
with a 10% iron content, as presented bv IChihara et al. I (l200l . 
Unfortunately, the 14jum complex is the only wavelength region 
where the enstatite features are not blended with forsterite emis- 
sion. This makes it impossible to study the enstatite iron content 
using other complexes. 

In the Galactic sample the 16 fim complex is clearly visi- 
ble in most sources, the outliers being EPLyr and IRAS 10174. 
The feature seen in IRAS 15556 is strongly deviating from the 
mean complex, and shows a stronger contribution of enstatite 
emission. IRAS 06338 shows strong emission of CO2 gas in the 
16fj.m region. 



As was already discussed in lGielen et alJ (l2008h . the strength 
of the I6fj.m feature seems to correspond to the emission of 
forsterite, shifted bluewards in central wavelength. Our new 
spectra follow this trend. This shift of the 16yum feature is also 
seen in the infrared spectra of protoplanetary discs (Juh asz et al.l 
2010), and is probably an effect of the adopted synthetic spec- 
trum of forsterite. 
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Fig. 6. Left: The lOyum complex for the Galactic sources, continuum subtracted and normalised. Overplotted in red the mean 
spectrum. The mass absorption coefficients of amorphous olivine and forsterite are plotted in blue and green. Right: Same as on the 
left, but for the LMC sources. The top panel shows the comparison between the calculated means for the LMC and Galactic sources. 
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Fig. 7. Comparison between the normalised and continuum- 
subtracted spectrum of one of our sample stars ST Pup and the 
laboratory spectra of ortho-enstatite (dot-dashed line) and clino- 
en statite (dashed line) with a 10% iron content, as presented 
bv lChihara et al.l (l2002h . The strong observed feature in ST Pup 
around 16//m is due to forsterite. 



4.3. The 19fim and 23fim complexes (17-21 jum and 
21 - 26 fim) 

This region shows two strong emission complexes, around 19 
and 23 yum , and a good agreement between the two samples is 
found (see Figs IA.2I and IA.3l l. Although the LMC sample has 
strong noise, the mean spectrum is very similar to the Galactic 
mean. 

The 19;um feature seems to be more pronounced in the 
Galactic sources, which could mean that the LMC sample 



is less crystalline, since the feature is mainly formed by 
forsterite emission. However, the strong noise level of the 
LMC sources could also hamper the detection of the feature. 
The 23 ijm feature is dominated by the emission of forsterite, 
and is clearly seen in most Galactic and several LMC stars. 
The mean Galactic and LMC complexes are again very sim- 
ilar. In the Galactic sources, the same outliers appear again: 
EP Lyr and HP 5296 1 , which have a very particular mineralogy 
(iGielen et al.l l2009b'): IRAS 10174, which is almost completely 
amorphous, and IRAS 15556 which shows no emission at 23 /zm. 

4.4. The 27 /jm and 33 complexes (25.5 - 30 fim and 
32 - 36 i^m) 

The two samples show a similar observed mean features, peak- 
ing around 27 and 33 ;um (see Figs. lA.4] and lA.5l ). 

The 27 yum complex peaks at the forsterite 27.5 fim feature, 
but with an additional shoulder around 29.3 yum, which is due to 
enstatite emission. In the Galactic sample, strong deviation can 
again be seen for IRAS 06338, which is clearly a source with 
atypical dust emission features. Some sources show a somewhat 
broader 21 fim feature, such as IRAS 15556, which could point 
to a larger enstatite contribution. The LMC sample is again com- 
promised by the strong noise, but the stars with strong observed 
emission feature do show a similar feature as observed in the 
Galactic sample. Only MSX 949 seems to deviate from the ob- 
served mean complex, with a very broad feature which peaks at 
29.5 Aim. 

Both samples show a clear 33 fj.m feature, due to the emis- 
sion of forsterite crystals at lower temperature. For the LMC 
sample the spectrum around the 33 fim complex is strongly ham- 
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Table 3. Overview of the adopted dust species. For each component we Hst its chemical composition, whether it has an amorphous 
(A) or crystalline (C) structure, density, adopted grain shape and grain sizes, and reference to the refractory indices used. 



Dust species 


Composition 


Structure 


Density 


Shape 


Grain size 


Reference 


Olivine 


Mg2Si04 


A 


3.71 g/cm' 


GRF 


0.1 


-2 


- 4yum 


Dorschner et al. (1995) 


Pyroxene 


MgSiOj 


A 


3.20 g/cm^ 


GRF 


0.1 


-2 


- 4yum 


Dorschneret al. (1995") 


Olivine 


MgFeSi04 


A 


3.71 g/cm^ 


GRF 


0.1 


-2 


- 4/jm 


Dorschner et al. (1995) 


Pyroxene 


MgFeSi206 


A 


3.20 g/cm^ 


GRF 


0.1 


-2 


- 4/jm 


Dorschner et al. (1995) 


Forsterite 


Mg2Si04 


C 


3.33 g/cm^ 


GRF 


0.1 


-2 


- 4/jm 


Servoin & Pirou (1973) 


Ortho-Ens tatite 


MgSiOs 


C 


2.80 g/cm^ 


GRF 


0.1 


-2 


- 4/jm 


Jaeger et al. (1998) 


Silica 


SiOj 


A 


2.20 g/cm^ 


GRF 


0.1 


-2 


- 4yum 


Henning & Mutschke (1997) 



pered by high noise, but the feature is still visible in the mean 
spectrum. A few sources have very strong emission at 33 yum, 
such as HV 12631, J044458, MACHO 78.6698.38, MSX949, 
SAGE 0543 10, and SAGE 050830. 



5. Full spectral model 

To study the characteristics of the silicate emission observed in 
these sources, we constructed a basic model to fit the full Spitzer 
wavelength range. The observed emission features will depend 
on the chemical c omposition o f the d ust, the grain sizes and 
the grain shapes. In lGielen et al.l (1200 8') we constructed a model 
that takes all the above properties into account. Note that a bug 
was pr esent in the modelling r outine used in Gielen et al.l (l2008[ 
l2009al) . which we describe in iGielen et al.i (2010) . For this pa- 
per, we also extended the routine to include an additional dust 
species, namely amorphous silica (SiOa). 

Assuming that the dust features are formed in an optically 
thin upper part of the disc, the spectrum can be approximated 
as a linear combination of dust absorption profiles. The model 
emission is then given by 

F, ~ aiKi) X (J^/3jB,(Tj)) + F,on, 
• j 

where is the mass absorption coeflncient of dust component / 
and a, gives the fraction of that dust component, B,\{ Tj) denotes 
the Planck function at temperature Tj and fij a scaling factor for 
the Planck functions. A sum of two Planck function s is also used 
to re present the continuum flux Fcom ■ Following iGielen et alJ 
(1200 8). we use two different dust and continuum temperatures, 
ranging from lOOK to lOOOK. 

The dust species we included are amorphous 
olivine/pyroxene (Mg2^Fe2(i-.v)Si04/Mg.,Fei_j:Si03), crys- 
talline olivine/pyroxene (forsterite/enstatite) and amorphous 
silica. Silica has different polymorphs, such as quartz, cristo- 
balite and tridim ite, with similar emission profiles (e.g. 
ISargent et al.ll2009 l). and we cannot rule out that some of these 
other polymorphs contribute to the silica fraction. To keep the 
number of free parameters to a minimum, we opted to use only 
amorphous silica in our modelling. In section 14.21 we showed 
that (part of) the enstatite content in our discs might be in the 
form of clino-enstatite with a 10% iron content. Unfortunately, 
the laboratory data of this enstatite species does not allow to 
calculate synthetic spectra for different grain sizes, so we opted 
to use the more commonly used iron-free ortho-enstatite. As 
discussed above (Sect. |4|, our study of the different complexes 
show that these dust species are present, and that there is no 
strong evidence for the presence of other dust species. 

Mass absorption coefficients for the different dust species 
are calculated from refractory indices in gaussian random fields 



(GRF) dust approximation (* Shkuratov & Grvnkol2005l) . The de- 
tails of the different refractory indices that we used can be found 
in Table [3] From our previous spectral studies we know that the 
observed emission features are reproduced using a non-spherical 
grain shape. Even though t he continuous distr ibution of ellip- 
soids approximation (CDE, Bohren et alj|1983l) is widely used, 
it is unfortunately only valid in the Rayleigh limit, and does not 
allow us to study grain growth effects. For this reason we prefer 
the GRF approximation. We also tested the distrib ution of hol- 
low spheres approximation (DHS. lMin et alJl2005h . but this did 
not result in a better fit to the observed emission features. 

To study the grain size distribution inferred from the mod- 
elling, we use three discrete dust grain sizes in the model: 0.1, 
2.0 and 4.0 //m. The emission features of grains with larger sizes 
be come too we a k to d istinguish from the continuum emission. 
In IGielen et alJ (l2008l) we already found that the presence of 
Mg-rich amorphous grains cannot be ruled out, and thus here 
also we use both purely Mg-rich amorphous silicates (x = 1) 
and amorphous silicates with an equal amount of Mg and Fe 
{x - 0.5). The ratio of magnesium and iron in the amorphous 
silicates mainly change s the peak positio n of the 10 and 18//m 
emission features ( Dorschner et al J lI995h . 

The best model was calculated using standard minimali- 
sation. Errors on the model parameters were calculated using a 
100 step Monte Carlo simulation with gaussian noise distribu- 
tion. Even though this model is only a first approximation, the 
model clearly succeeds in giving an overall good fit to the ob- 
served spectra (see Figs. lA.6llA.10l l. 

For 6 sample (Galactic) sources the Spitzer spectrum only 
starts at 9.9 /vm, which means we lack information on the dust 
composition in the 10 /im wavelength range. Since this could in- 
fluence the derived dust parameters, we depict these sources in a 
different color in our correlation plots (Figs[8l lA.12l ). 

5.1. Results 

We find that, for most sources, the dust is dominated by large 
grains. We define the mass-weighted mean grain size of the dust 
as 

with fl, the grain size, and m, the mass fraction of dust in that 
grain size. For 47/57 of our sample stars, the mean grain size is 
larger than 2 yum (see Fig. |8]l. 

Since our model routine uses three different grain sizes, we 
can use them to determine a grain size distribution. The Spitzer 
spectra probably only trace the upper layers of the disc, and so 
the calculated distribution could not be valid for the entire disc. 
The grain size distribution is usually approximated by n{a) oc 
a'\ with n{a) the number of grains with grain size a, and p a 
power-law ind ex. For the interste llar medium a value of p - 
-3.5 is found (iMathis et al.lll977h for typical ISM grains up to 
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0.3 ^m , rolling over exponentially for larger grains (IZubko et alJ 
12004 V To calculate the number of grains in a given grain size, 
we compute the mass fraction of these grains from our modelling 
and divide it by the corresponding volume of the grains. We then 
normalise all the grain numbers, such that n(0.1 fim) = 1. The 
results of this calculation can be seen in Figure|9] We find a good 
fit to our results is achieved with a power-law index p - -I .Sg j, 
for grain sizes between 0. 1 and 4 fim . It is clear that our grain 
size distribution is not ISM like, larger grains are much more 
abundant. 




12 3 4 

Meon Grainsize [jLtm] 

Fig. 8. The mass fraction in crystalline grains versus the mean 
grain size of our spectral modelling. Galactic sources are given 
in red plus signs and LMC sources in blue diamonds. The ma- 
genta symbols depict Galactic sources for which the infrared 
spectra only start from 9.9 fim. The numbers correspond to num- 
bers given in Tables[T]and|2] 



From Figure IA.12I it is clear that for about half the stars 
the crystalline grains are larger than the amorphous grains. For 
nearly all stars the mean grain size of the crystalline grains lies 
above 2 yum, whereas the amorphous grains show a larger spread 
in grain sizes. This is in contrast to what is found for the dust 
in discs around Herbig Ae stars, where the cryst alUne grains 
are s ignificantly smaller than the amorphous grains ( Juhasz et al.l 
1201 0). We do not find any correlation between the size of crys- 
talline and amorphous material. It is unclear what causes this dif- 
ference in grain size between the amorphous and crystalline dust. 
An effect that could come into pla y here is the appa rent spectral 
signature of large dust aggregates. |MineLa n (l2008h showed that 
aggregates with a very low abundance appear spectroscopically 
as very small grains, while more abundant materials appear spec- 
troscopically to reside in larger grains. Since for our sources the 
amorphous dust is in most cases more abundant than the crys- 
talline dust, this could mean that the amorphous grains reside in 
large fluffy aggregates, which have spectral signatures that are 
very similar to those of small g rains (iMin et ail2606li2008h . 

Similar to what is found in I Juhasz et al.l (1201 Ol) we find that 
the size of the enstatite grains is on average slightly larger than 
that of the forsterite grains. There seems to be a weak trend 
between the crystallinity and the mean size of the crystalline 
grains: sources with a higher crystallinity have on average larger 
crystalline grain sizes (see Fig. lA.12T i. 
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Fig. 9. The normalised number of grains versus the adopted grain 
sizes. The triangles represent all the sample stars. The solid line 
gives the best power-law distribution to the mean of all stars. 
The dashed lines represent different power-law indices, given for 
comparison. 



For discs around young stars a strong correlation is found be- 
tween the mean grain size of the amorphous grains and the disc 
flaring. This disc flaring is determined by the ratio of the 24 //m 
and 8/im flux. Juhasz et al . (2010) find that sources with flat- 
ter discs have larger amorphous grains in their disc atmosphere. 
This trend is not seen in our sample. Our sample sources seem to 
be centred around ^24/^8 - 0.83, which shows that these discs 
are not strongly flared and that there is no large spread in disc 
flaring. Of the 57 sources, 8 show higher values of the disc flar- 
ing, going from F2aIF% - 2 up to 4.3. We find no correlation 
between the disc flaring and any other dust parameter. 

Figure [8] also shows the high crystallinity fraction derived 
from the Spitzer modelling. The crystallinity can reach values of 
60%, which is among the highest seen in astronomical environ- 
ments. High values of crystalline dust ar e also found for proto- 
planetary discs around young stars (e.g. Bou wman et al.lB008l: 
ljuhasz et al. 2010), where crystalline fractions up to about 30% 
are found. For the crystalline dust, forsterite is almost always the 
dominant species: the forsterite fraction of the crystalline mate- 
rial has values between 20 and 100% (Fig. lA.lll i. There does 
not seem to be a strong correlation between the crystallinity and 
the forsterite/enstatite fraction of the crystalline material. The 
same holds for the forsterite/enstatite fraction of the crystalline 
material and the mean grain size of the crystalline grains (see 

Fig. EH). 

The derived silica fractions are of the order of 5%, but can go 
up to 20%. We find no correlation between the silica mass frac- 
tion and any other dust parameters. Thermal annealing of amor- 
phous dust produces both forsterite and silica, and so a relation 
is to be expected if this process is responsible for the crystallisa- 
tion. 

In our modelling we used amorphous silicates with an equal 
Mg-Fe content and pure Mg-rich grains. We find that 25/57 
sources show a dominance of iron-free amorphous dust, whereas 
the other sources are clearly dominated by Mg-Fe amorphous 
silicates. On average we find that the Galactic sources have a 
slightly higher fraction (53%) of the Mg-Fe rich amorphous sil- 
icates, whereas the LMC sources have a higher fraction (56%) 
of purely Mg-rich amorphous silicates. The derived differences 
are very minimal, so we cannot make strong statements on the 
iron content of the amorphous silicates. For the sources where 
we find a high fraction of iron-free dust, the iron grains could be 
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Stored as metallic inclusions in the grains, which would be very 
hard to detect. 

We do not find any correlation between the derived dust pa- 
rameters and central binary parameters such as the effective tem- 
perature or the orbit. 

5.2. Atypical sources: intruders? 

Some sources clearly deviate from the mean observed spectrum, 
by showing no crystalline grains (IRAS 10174 and J05133) 
or carbonaceous molecules and/or gas emission (EPLyr, 
HD 52961, IRAS 06338 and IRAS 13258). Other sources, such 
as MSX949, show less obvious differences, but are still not 
reproduced as well by the model as the other sources. Since 
these sources were mainly selected on the basis of their infrared 
colours, we cannot exclude that non-post-AGB disc sources 
are present in the sample. Possible intruders could be young 
stars with protoplanetary discs, red super giants or AGB stars. 
However, for most sources we have additional observations of 
the central star, such as optical spectroscopy, which corroborate 
their post-AGB evolutionary phase. We discuss several doubtful 
(or anomalous) cases below. 

5.2.1. SAGE 050830 

For SAGE 050830, the optical spectra show some evidence for a 
carbon-rich chemistry ( van Aarle et al. 201 1). However, the pho- 
tometry and infrared spectral information for this source is sus- 
pected to be contaminated by a foreground star of spectral type 
AO-IIV. Unfortunately, the angular resolution does not allow us 
to discriminate between the A star and the carbon star as the 
identification of the Spitzer source. Still, if the carbon-rich spec- 
trum truly belongs to the Spitzer source, the strong oxygen-rich 
spectrum seems surprising. This source is one of the more crys- 
talline objects of our sample, and even has the most extreme 
lOyum feature-to-continuum ratio of all sources! 

The carbon-rich classification of the central star, together 
with the presence of crystalline silicates in its circumstellar envi- 
ronment would make this star an ideal candidate to be a silicate 
J-type carbon star. These are carbon-rich AGB stars, but with a 
very low '^C/'^C ratio and detection of crystalline silicates i n 
their infrared specti-um (iLlovd Evanslll990tlAbia & Isernll2000h . 
The sources are believed to be binary stars, w ith an unseen 
companion, surrounded by a circumbinary disc (Morris 
[jura & Kahane 1999; Yamamura et al. 2000; Deroo et al. 2007). 
This scenario could explain the dual chemistry, since the disc 
could then be formed while the central star was still oxygen rich, 
and has now evolved to be carbon rich. However, it does not ex- 
plain the low '^C/'^'C ratio, usually seen in J-type silicate carbon 
stars. 

Unfortunately, due to the confusion with the foreground star, 
we cannot determine the stellar parameters, which would shed 
light on the evolutionary status of this object. Also, the low- 
resolution optical spectrum does not allow to determine the 
'^C/'^C ratio of the carbon star, hence corroborating the J-type 
nature. 

5.2.2. J 051 33 and IRAS 101 74 

Two of our sample sources, IRAS 10174 and J 05133, clearly 
deviate from the rest of the sample by showing no strong evi- 
dence for crystalline features in their spectra. However, our mod- 
elling shows that purely amorphous silicate dust is not sufficient 



in reproducing the observed features, especially around 13 jum. 
The spectra of these stars are actually very similar to the ob- 
served spectra of AGB outflows, characterised by small amor- 
phous grains. In these sources, emission from additional dust 
species, such as alumina (AI2O3), can influence the 13/vm re- 
gion. To see if alumina could also be present in these sources, 
we remodelled the spectra, now including alumina grains. For 
both stars we find an improvement when including AI2O3, es- 
pecially for J 05133. For IRAS 10174 the improvement is only 
minor, with an amount of alumina in the new model of 2%. 
However, for J 05 1 33 the fit is improved drastically when includ- 
ing 30% alumina (Fig. [TOll. The bulk of the other dust (-50%) is 
stored in small 0. 1 fim Mg-rich olivine in this new model. Less 
then 10% of the mass fraction of dust is in crystalline form, and 
then mainly forsterite. This type of dust composition is more in- 
dicative of an outflow and not a disc. 

The optical spectrum of J 05133 dvan Aarle et al.l[2"01 1[) also 
indicates the peculiar nature of this source. The spectrum points 
to a F8-G0Ip spectral classification, but shows very strong HI 
(6563 A) and He I (5876 A) emission, and broad Ca II absorption 
lines, which are not expected in a star of this type, but point to 
the presence of a hotter source. One possibility would be that the 
system is actually a binary with an unseen hot companion. The 
SED modelling gives, for this source, a luminosity of around 
17 000 Lq. This, combined with the spectral type as derived from 
the optical spectrum, shows that the source cannot be an AGB 
star, but also shows that it is probably not a post-AGB disc source 
as normally understood. 



J051333. 74-663419.1 
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Fig. 10. The results of our modelling of J 05133.74, with the ad- 
dition of alumina grains. The observed spectrum (black curve) 
is plotted together with the best model fit (red curve) and the 
continuum (black solid line). Crystalline silicates are plotted in 
green, amorphous sihcates in magenta and alumina in blue. 



The unusual chemistry of IRAS 10174 compared to the rest 
of the sample made us re-investigate the optical spectra for 
this source. The source was or i ginally classified as a post- 
AGB source dLlovd EvansI ll999l:lDe"Ruvter et al.1 |2006ft from 
a blue spectrum which suffered badly from dust extinction 
in the violet. A new optical spectrum (6400 - 9000A), taken 
with the Cassegrain grating spectrograph of the 1.9 m RadclifFe 
Telescope at the South African Astronomical Observatory, 
shows that the star is a luminous supergiant, of type G8Ia- 
O. This is in better agreement with the observed chemistry of 
the infrared spectrum. Note that even though the star is not 
a post-AGB disc source, the binarity of the star is confirmed 
(iMaas et al ] l2003h . 
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Even though some non-post-AGB disc sources might be 
present in the sample, these sources do not change the overall 
conclusions of this study. The removal of these sources does not 
introduce correlations between the different stellar and dust pa- 
rameters, that are currently not observed. 

6. Discussion 

For all sources, the infrared spectra are dominated by emis- 
sion features due to oxygen-rich dust species. For some 
of the LMC sources (such as J051159.ll, J053605.56, and 
SAGE054310) this is surprising, since they have luminosities 
which would put them in the peak of the carbon star lumi- 
nosity functio n (Stancliffe et al. 2005: lGroenewegen et al.ll2007t 
ISrinivasan et al.li201 li) . Optical spectroscopy for these sources 
does not point to a c arbon-rich chemistry of the central star 
(Ivan Aarle et al.ll20ll1) . This shows that the AGB evolution for 
these binary post-AGB stars was shortcut, possibly under the in- 
fluence of strong binary interaction, preventing them to evolve 
into carbon stars. 

Our study shows that even a relatively simple model suc- 
ceeds in reproducing the observed infrared spectra. The model 
assumes only two dust temperatures, and uses the same dust 
abundances for the cool and warm dust. This shows that the dust 
in the disc is relatively well mixed, in that cooler and hotter re- 
gions in the disc have a similar dust composition. This is very 
different from the results for protoplanetary discs around young 
stars, where a difference in dust composition is needed for the 
inner and outer disc regions ( Juhasz et al. 2010). 

The strong observed crystalline bands at longer wavelengths 
show that at least a significant fraction of the crystalline grains 
are located at cooler temperatures. For some sources our model 
even underestimates the forsterite flux at 33 //m, showing that, 
for some sources at least, the forsterite fraction between the two 
temperatures might not be evenly distributed, but dominated by 
the cooler temperature. The problem of reproducing the features 
at longer wavelengths might also be due to the adop ted synthetic 
spectra of forsterite. As can be seen in Fig. 3.14 of Gielen et al.l 
(l2008h . the 33 yum is best reproduced by DHS grain shapes. The 
GRF grain shape gives a feature which is much broader and flat 
topped. However, since on average our features were slightly 
better reproduced with GRF shapes, we used this approximation 
in our modelling. Another effect that can influence the observed 
features in the optical depth. At different wavelengths we would 
look at different depths in the disc, with a different temperature 
distribution. At longer wavelengths, we would then look deeper 
in the disc, where the cooler temperature might enhance the fea- 
tures at these longer wavelengths. To study this effect in detail, 
a full radiative transfer model is needed, which goes beyond the 
scope of this paper 

Since the exact formation mechanism of these circumbi- 
nary discs is still uncertain, it is difficult the relate the differ- 
ent observed dust characteristics to disc evolution. One possibil- 
ity is that the discs are formed after a common-envelope phase, 
with some dust formation already forming in the outflow phase. 
Crystalline grains can then be formed directly out of the gas 
phase, at high temperatures (Gail 2004; Petaev & Wood 2005). 
But this will probably not give rise to the very high amounts of 
crystalline material we see, which shows that another crystallisa- 
tion process is still a ctive afterwards, such as thermal annealing 
(IWooden et al.ll20d 5'). Another formation mechanism is Roche- 
Lobe overflow through an outer Langrangian point, where the 
material is already confined to the midplane. Here one could 
expect the dust at the hot and denser inner regions to be more 



efficient in producing crystalline species. Both gas-phase con- 
densation and thermal annealing might be important to explain 
the high crystallinity in these discs. 

The condensation models predict forsterite to condense first, 
followed by the formation of enstatite through reactions between 
forsterite and SiOa gas. In contrast to what is found for proto- 
planetary discs, our results show that forsterite is almost always 
the dominant crystalline dust species. This could point to a de- 
viation from equilibrium conditions during condensation. If the 
forsterite grains reach large grain sizes quickly, the formation of 
enstatite might be complicated, since it will become increasingly 
harder to infuse SiOa in the forserite lattice. The resulting dust 
may the be in the form of a large forsterite grain, surrounded 
by a small layer of enstatite. The formation of enstatite can be 
further weakened if the material is allowed to cool very quickly 
after the condensation of forsterite. The formation of forsterite 
through annealing is especially efficient if the starting material 
has an olivine stoichiometry. A high forsterite fraction would 
then go together with a higher pyroxene fraction of the amor- 
phous material, which is not supported by our results. 

Sinc e crystallisat ion requires high temperatures above 
lOOOK dFabian et al.l 12000). one would expect the crystalline 
dust to be confined to the hot, inner regions of the disc. This is 
in clear contrast to our findings of cool crystalline material and a 
homogeneous dust composition throughout the disc. This shows 
that mixing must be efficient is transporting the crystalline ma- 
terial to cooler regions which were initially dominated by amor- 
phous material, or a crystallisation process at lower temperatures 
is occurring in the discs. Molster et al. (2002a b) already showed 
that the crystallinity fraction in disc sources is much higher than 
that observed in typical outflow sources. This shows that the 
crystalline component in the disc sources is most likely deter- 
mined by subsequent dust grain processing in the discs, and not 
by cooling processes in the outflow of the material forming the 
discs. 

Dust formation models also show that iron will preferably 
condense out as metallic iron, rather than be included in silicate 
formation. This could explain the presence of Mg-rich amor- 
phous silicates in our results. The grain sizes of the crystalline 
dust, formed through condensation, will not be correlated with 
the grain sizes of the amorphous material, which is in line with 
our results. If the crystallisation occurs through annealing, we 
would expect a relation between the initial amorphous material 
and final crystalline grains. This does not explain the observed 
difference in crystalline and amorphous grain sizes, unless a sub- 
sequent process can be invoked that would grow the crystalline 
material, but not the amorphous dust. Our results on the differ- 
ence in crystalline and amorphous grain sizes is again in clear 
contrast to what is found for the dust in protoplanetary discs, 
where the crystalline grain s are found to be sig nificantly smaller 
than the amorphous grains (l Juhasz et alj20l d). Clearly, different 
dust processes are responsible for the grain growth and crystalli- 
sation in the discs around young and evolved stars. 

The derived large grain sizes show that there seems to be an 
efficient removal of the smallest grains. The question remains 
whether this lack of small grains is an effect of grain growth 
(Dullemond & Dominik 2004) or whether the initial grain pop- 
ulation already consisted of large grains. In that case the small- 
grain fraction could be a result of grain collision and subsequent 
break-up. An effect which might also be important to the ob- 
served grain sizes is the strong radiation of the central source. 
The central post-AGB stars are highly luminous, and radiation 
pressure could be responsible for the removal of the smallest 
grains in the upper layers of the disc. Since our results show 
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that the amorphous grains tend to be smaller than the crystalline 
grains, radiation pressure might be (partly) responsible for the 
large fraction of crystalline grains observed in the upper layers 
of the disc. 

Surprisingly, we find no correlations between the derived 
dust parameters, such as crystallinity, grain size and abundances. 
Also, no correlation between the dust parameters and parameters 
of the central binary system is found. The lack of correlation 
raises the question whether the optically thin upper layers traced 
by the Spitzer spectra are a good representative of the global dust 
composition. 

Except for the amorphous silicate dust, we find no evi- 
dence for the presence of dust species usually associated with 
AGB outflows or single-star post-AGB shells, such as simple 
oxides or Al/Ca-bearing dust species. The theoretical oxygen- 
rich dust condensation sequence for dusty outflows starts with 
the formation of alumina (AI2O3) around 1760 K, followed by 
formati on of gehleni t e (CaiAliSiO?) at s lightly lower temper- 
atures (lTielenslll990l: iTielens et al.lll998h . Further interactions 
with magnesium will produce species like spinel (MgAl204), 
akermanite (Ca2MgSi207), diopside (CaMgSi206) and finally 
anorthite (Ca2Al2Si208) around 1360 K. A second condensation 
sequence, involving mostly magnesium and silicon, starts with 
the formation of forsterite around 1500 K, followed by enstatite 
around 1300K. Only at temperatures below the glass temper- 
ature can iron interact to form amorphous iron-containing sili- 
cates. Not only the temperature plays a role, also the densities 
involved will determine which dust species can be formed. 

From observations it is found that AGB stars start by form- 
ing Al- and Mg-rich oxides in their outflows, followed by an 
increase of amorphous silicates bands with incre asing mass-loss 
rate, which start to grow on the Al-rich oxides dLebzelter et alj 
I2OO6I). A similar tren d is seen in outflows of red supergiants 
dVerhoelst et alJ2009h . A similar scenario might explain the lack 
of Al/Na/Ca-rich dust species in the discs around the post-AGB 
stars. Since the photosph eres of the central stars are strongly de- 
pleted in these elem ents dMaas et alJl2005t iHrivnak et aljboOSt 
iGielen et alj2009bh . we know these refractory elements must be 
present in the disc. Of course, since densities associated with 
these discs are much higher than for typical outflows, and dust 
might be subject to a different temperature gradient, which could 
result in a different condensation sequence to that observed in 
AGB stars. 

7. Conclusions 

We analysed the Spitzer infrared spectra of 33 Galactic and 24 
LMC (candidate) post-AGB binaries surrounded by a dusty cir- 
cumbinary disc. For nearly all Galactic sources, previous studies 
have already confirmed the binarity and post-AGB status. The 
LMC sources were taken from a list of probable post-AGB disc 
candidates. Our main focus was to determine the dust composi- 
tion of the discs, but also to look for possible differences between 
the Galactic and LMC sample. Our study shows that: 

- The Spitzer spectra are all dominated by emission features of 
oxygen-rich dust species, namely amorphous and crystalline 
silicates of olivine and pyroxene stoichiometry. 

- The observed silicate dust has a high crystallinity factor: 
most sources have crystalline mass fractions between 20 - 
60%. 

- Most of the dust is stored in larger grains (> 2 yum ). This 
results in an average grain size distribution of n{a) <x a^^\.2 
for grain sizes between 0. 1 and 4 //m . 



- We find no correlations between the dust, stellar, and/or or- 
bital parameters, which makes it difficult to constrain the 
dust grain processes that are causing the observed dust prop- 
erties such as the grain sizes and crystallinity. 

- We find no differences between the dust parameters of the 
Galactic and LMC sources. 

- Although the observed spectra are very similar to those 
of protoplanetary dies, we find evidence for a fundamen- 
tal difference in the dust processing occurring in the two 
disc types, more specifically in the homogeneous dust com- 
postion throughout the disc, the observed degree of crys- 
tallinity, the crystalline grain sizes and, the strong dominance 
of forsterite in the crystalline grain fraction. 
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Appendix A: Figures and Tables 
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Fig.A.l. Left: The 14//m complex for the Galactic sources, continuum subtracted and normalised. Overplotted in red the mean 
spectrum. The mass absorption coefficients of forsterite and enstatite are plotted in green and blue. The high noise level on the LMC 
sources did not allow for a mean spectrum determination. Right: Same as on the left, but for the 16/im complex. 




Fig.A.2. Left: The 19 fim complex for the Galactic sources, continuum subtracted and normalised. Overplotted in red the mean 
spectrum. The mass absorption coefficients of forsterite and enstatite are plotted in green and blue. Right: Same as on the left, but 
for the LMC sources. The top panel shows the comparison between the calculated mean for the LMC and Galactic sources. 
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Fig. A.4. Same as Figure IAT2I but for the 27 fim complex. 



extragalactic post-AGB disks., Online Material p 2 



Goioctic 




JK 




r MSX949 


^ SAGE05431 




m l4/2B62 




TO [jApH082.840Syi 5 


in 


35274 






«|\ J0^./.3.1S 


^ MACH07S-6698-3S 












ki 




W HV2281 






i/V 

V SAGE05083d / 


n HV12631 1 Vl 




if 522 20. 8 7 






□5 RBI 


21 23 25 27 21 23 25 27 2 1 23 25 27 21 
Wavelength [^tm] 


33 


5 2 


Vlmc 

Goiactic ^ 












HV5829 ^ 




MSX949 ^ 


SAGE054310 ' 


HV2862 ^ 






UsA^^052747 




1^2522 1/^ 


JDqqi43.i8 ^ 




m 








SAGE051453 ^ 


HV2281 ' 




J053605.56 ^ 


SAGE050830 ^ 


/\ 

1'hvi2B31 ^ 
.5 27.5 29.5 25 


k% 

.5 27.5 29.5 25 


J05222D-B7 ^ 
.5 27.5 29.5 25 


(nGCI aaSSAGEIRSl 
.5 27.5 29.5 



C. Gielen et al. 



Silicate features in Galactic and extragalactic post-AGB 



disks., Online Material p 3 




32 33 34 35 36 32 33 34 35 36 32 33 34 35 36 32 33 34 35 36 32 33 34 35 36 32 33 34 35 36 32 33 34 35 36 32 33 34 35 36 



Fig. A.5. Same as Figure lA!2l but for the 33 jim complex. 
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Fig. A.6. Best model fits for our Galactic sample stars, showing the contribution of the different dust species. Top: The observed 
spectrum (black curve) is plotted together with the best model fit (red curve) and the continuum (black solid line). Forsterite is 
plotted in green, enstatite in blue, silica in cyan and amorphous olivine and pyroxene in magenta. Small grains (0. 1 fim) are plotted 
as dashed lines and larger grains (2 and 4 yum) as dotted lines. Bottom: The normalised residuals after subtraction of our best model 
of the observed spectra. 
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Fig. A.7. Same as Fig. lA.6l 
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Fig. A.8. Same as Fig. lA.6 
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Fig.A.9. Best model fits for our LMC sample stars, showing the contribution of the different dust species. Top: The observed 
spectrum (black curve) is plotted together with the best model fit (red curve) and the continuum (black solid line). Forsterite is 
plotted in green, enstatite in blue, silica in cyan and amorphous olivine and pyroxene in magenta. Small grains (0. 1 jim) are plotted 
as dashed lines and larger grains (2 and 4 /im) as dotted lines. Bottom; The normalised residuals after subtraction of our best model 
of the observed spectra. 
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Fig. A.IO. Same as Fig.lAigl 
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Fig. A.ll. Left: The mass fraction in crystalline grains versus the forsterite fraction in the crystalline material. Right: The mean size 
of crystalline grains versus the enstatite fraction in the crystalline material. Galactic sources are given in red plus signs and LMC 
sources in blue diamonds. The magenta symbols depict Galactic sources for which the infrared spectra only start from 9.9 fim. The 
numbers correspond to numbers given in Tables [1] and |2] 
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Fig. A.12. Left: The mean size of crystalline grains versus the mean size of amorphous grains. Right: The mean size of the crystalline 
grains versus crystallinity fraction. Galactic sources are given in red plus signs and LMC sources in blue diamonds. The magenta 
symbols depict Galactic sources for which the infrared spectra only start from 9.9 /vm. The numbers coiTespond to numbers given 
in Tables [1] and m 




Fig. A.13. The mean size of the amorphous grains versus the disc flaring, determined by the /^24/^^8 flux ratio. Galactic sources are 
given in red plus signs and LMC sources in blue diamonds. The magenta symbols depict Galactic sources for which the infrared 
spectra only start from 9.9 yum. The numbers correspond to numbers given in Tables[T]and|2l 
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Table A.l. Best fit parameters deduced from our full spectral fitting. Listed the dust and continuum temperatures and their relative 
fractions. 



N° Name 



Tdustl 

(K) 



1(KF 
200 

100, 

200.; 

426 
200.f" 
200.°- 
179 
200.°- 
209.f 
214f 

wo.i 

111 
100.; 
100.; 
200.0- 

218-!^ 
320 
200.°- 
100.0- 
136 
201 
100 
200, 
100 
277, 
101 

203.''9- 
213 



Tdustl 

(K) 



200":^ 

lOOO.i!- 
1000.^- 
704.i" 

771 204. 
' ' ^-153. 

680.Jg; 
1000.0- 
600.°- ■ 

627.2f- 

699 

300.jJ 
426.^5. 

200.°- ■ 
200.°- 
500.°- 
952.^^- 

137. 
64. 

500.« 

200.°- 

240."- 
414^81. 

200.°- 
799.; 

200.!; 

576 
998 



Fraction 



90" 

90° "° 

■ !!-lB 
90° °° 

■ !!Q!! 
90°-°° ■ 

70°-^° ■ 

"•'"0.30 

90° °° 
90°-°° ■ 

60°^° 

■ !!-5Q 

90° °° 
"■^"0.00 

0-90°;j;°- 
0.80°i°. 
0-90°;S!! 

0.80°:'° . 

0.80; - 

20° '° - 

9.19. 

90° °° - 
"■^"0.00 

0.90°;?°- 

0.80°;1°. 

90° °° ■ 

■ !!Q!! 

90°-°° ■ 
Q-IQ 

0.80°- ° ■ 



-0 10™" 
"•^"0.00 

- 10°-'° 

Q-!!Q 

- 10°-°° 
Q-flQ 

■ 10°-°° 
Q.OQ 

■ 30°-"° 
"•-".'0.20 

- 10°-°° 
^-'•"'0.00 

- 10°-°° 
^-'•^"0.00 

- 40°-^° 
'■■'.^"0.30 

- 10°-°° 
"•^"0.00 

■ 10°-°° 
"•^"0.00 

- 20°-''° 
"•^"0.10 

- 10°-°° 
"•^"0.00 

■ 0.20°-?° 



TcontX 

(K) 



200.°- 

200. 
100. 

200.1, 

200.°- 
178.^f 
277.23- 
100.°; 
200.° 
376.24- 
100.°;" 
100.°- 
100. 
100. 
100.1, 
200.° 
200.°- 
200.°- 
100. 
100. 
200.S, 
500.° 

200.; 

304.'''- 

200.°- 

200.°- 

200.°- 

158.«- 

200.° 

100.°; 

300.° 

300.°- 

232.^?- 

203.''" 

304.ffl- 
200.° 
203. 
520.: 
431. 



Tcont2 

(K) 



996^ 
500.° 
800.°- 
1000.°- 
500.°- ■ 
500.°- 
564 ^f'^ 
859 
600.° 
700.°- 
976.2*- 
423.^«: 
603.''»- 
500.°- 
689.|J- 

500. ° 
513.'J- 
600.°- ■ 
600.°- 

600. °- 
464.|: 
700.° 
900.°- 

601. °- 
850.| 
596.f„ 

501. ? 
472.28- 
776 '3° 
617.?«- 
500.° ' 
800.°- 
900.°- 
605.Jf- 

504.!"" 
Tie. 

1%. 

■44. 



Fraction 



o.90g;r 
0-98 ; ■ 

99°-°° ■ 
Q-QQ 

93°-°° ■ 
940 00 . 
"■^^0.00 
94° °* ■ 

94°°' 
".'^^o.oi 

95° °° ■ 
93° °° 

"■^■^0.01 

91°°' 
"■^'0.01 

98° °' ■ 
".^°o.oo 

96°-°2 

Q.QQ 
QQO-OI' . 

QQfl 
0.99°°° 

95^-^''' ■ 

95° °' 

92° °° ■ 
"■'^0.00 

90° °° 

■ fi-Qfl 

97°-°° ■ 

993 

98° °° 

999 

96° °° 
"■^"0.00 

0.91°°° 

948:8* . 

99) 

0.91°-°' ■ 

098:84 . 

90° °° 
"■^"0.00 

94° °' 

"■^^0.04 

97° °' 
"■^'0.01 

93° °' 

95° °° 
".'^-^o.oo 

91°°° 
"■^'0.00 

97° °° ■ 
"■^'0.00 

93° °2 
"■'-'0.03 

O.o^o oo ■ 
0.81°;'° ■ 
0.87; - 
0.90; ■ 
92°-°' ■ 

".'-^0.54 

91°°' 
"."0.01 

95°°' ■ 

.8:8^ 

1.00 

93° "' ■ 

■ Q-Q? 
90° " 

97°-°° ■ 

99Q 
98° °' 

991 

"■^^0.02 
94°-°' - 
0.82f 
0.95; 
0.89°°' 

89^-^^ 
9H 
900.02 

90^-"^^ 
"■^"0.00 

0.88°»2 



-0.10°:° 

- 02°-°° 

Q.QQ 

■ 01°-°° 

QQO 

■ 07°-°' 

O.QQ 

- 06°-°° 

- 06°-°2 
^^•^^'■'0.04 

- 06°-°' 

- 05°-°' 

^-'•^-'.-'0.00 

- 07°-°' 

• 09°-°' 

• 02°-°° 

■ 04°-°° 

Q.Q2 

■ 01°-°° 

QQO 

■ 01°-°' 

QQO 

- 05°-°' 
■0.05^ 
■0.08JJ;°JJ 

- 10°-°° 

Q-QQ 

■ 03°-°° 

Q-QQ 

- 02°-°° 

Q-QQ 

- 04°-°° 

- 09°-°' 

Q-QQ 

- 06°-°° 

QQO 

■ 09°-°' 

Q.QJ 

■ 01°-°' 

"•"^0.00 

- 10° °° 
^-'•"-'o.oo 

- 06°-°* 

- 03°-°' 

- 07°-°' 

0-Q3 

- 05°-°° 
"."-'0.00 

- 09°-°° 
"."'^o.oo 

- 03°-°° 
"•"■^0.00 

■ 07°-°3 
"•"'0.02 

- 15°-°° 

".'-^O.OQ 
. O.I9O-2O 

- 13^-'^ 

9-9i 

- 10°-°° 

■ 0.08f 

■ 25" 

- 09"-°' 

- 05°-°2 

Q-Ql 

- 10°-°° 

QQO 

- 07°-°2 

QQl 

- 10°-°° 

Q-Q4 

- 03°-°° 

Q-QQ 

- 02°-°" 

- 0.06] 

- 0.06' 

- 05"-°° 
"•".^0.00 

- 11°-°" 
"."0.01 

- 12°-2* 

005 

- 10°-" 

Q-Q2 

- 10°-°° 

■ 12°-°3 
"•'^^0.02 

- 05°-°° 



1 EPLyr 

2 HD131356 

3 HD2 13985 

4 HD52961 

5 IRAS05208 

6 IRAS06034 

7 IRAS06072 

8 IRAS06338 

9 IRAS09060 

10 IRAS09144 

11 IRAS09538 

12 IRAS 10174 

13 IRAS 11000 

14 IRAS 13258 

15 IRAS 15556 

16 IRAS16230 

17 IRAS 17038 

18 IRAS 17233 

19 IRAS 17243 

20 IRAS 17530 

21 IRAS18123 

22 IRAS18158 

23 IRAS 19125 

24 IRAS 19 157 

25 IRAS20056 

26 RUCen 

27 SA0173329 

28 STPup 

29 SUGem 

30 SXCen 

31 TWCam 

32 UYAra 

33 UYCma 

34 HV 12631 

35 HV2281 

36 HV2444 

37 HV2522 

38 HV2862 

39 HV5829 

40 HV915 

41 J044458. 18-703522.8 

42 J045242.93-704737.4 

43 J050143. 18-694048.7 

44 J05 1159. 11-692532.8 

45 J051333.74-663419.1 

46 J052220.87-655551.6 

47 J053605.56-695802.9 

48 J054312.52-683356.9 

49 MACH078.6698.38 

50 MACHO82840515 

51 MSX949 

52 NGC1805SAGEIRS1 

53 SAGE050830 

54 SAGE051453 

55 SAGE052707 

56 SAGE052747 

57 SAGE054310 



•1. 

'•If- 
171.*f- 

206.^f 
919437. 

200.°- 
189.^1 
301.'f 
323 

nlf 

117. ^ 
250.S- 
230.^°5- 
209.f 

243.,'" 

100.; 

204.,'f- 
226.j' 

200.; 

183. ; , 
200.^-'' 
125.|- 
208.f- 
100.°; 
242.88. 

218.\ 

184. 'f 
203.2]"- 
248.^^; 
197 



'•127. 
lO. 

?0. 



563.. 



24. 

po. 



487 25- 

° ' -90. 

506.g«^ 

990.| 

400.° 

869.«; 

1000.° 

382.*°;- 

505.f^. 
916 

613. 

"^83.1^^; 
783.' „; 
452.^^ 
712 2"'- 

' -'^-104. 

1000.;;- 

709.1?°- 

990.1*- 

986.- 

812.11 
950. °3 
533 ^2": 

-'"^•74. 



200. 

582.2«; 
418.^*- 

«:| 



60° 3° 

90° °° 

■ Q-QQ 
90° °° 

■ Q-QQ 
0.90°-°° 



0.90°- 
0.90°- 



J:8§ 
"0.00 
0.80°1° 

0.80° 

0.80! 



0.10 

lO.lO 



0.70; 
0.70: 



'0.10 

lO.lO 



'0.30 

o.9o!!;S!! 

0.80-° 
90°'° 

".^"o.QQ 

50°-*° 

Q.3Q 

60°-^° 

9i9 
90° °° 

■ Q-AQ 

90° "° 

■ Q-IQ 
90° °° 

o.9o!!;S!! 

0.80°;'° 

00° "° 

■ Q-QQ 
90° °° 

■ Q-JQ 
10°^° 

■ Q-IQ 
90° °° 

■ Q-QQ 
90° °° 

■ Q-IQ 
90° °° 

■ Q-IQ 
90° °° 

Q.QQ 

90°-°° 

Q-lQ 

0.70i|;*« 

0.80!J;i!J 
0.80°;'° 

90° °° 

■ Q-IQ 
90° °° 

o.9o!i!!i 



465. 



0.90; 



8:^ 

'o.no 



0.20; 

0.80; 

10°-°° 
".'"0.00 

■ 10°-'° 
".'"0.00 

• 20°-'° 
".^"0.10 

- 10°-°° 
Q-OQ 

■ 10°-'° 
Q.QQ 

■ 20°-'° 
Q-iQ 

- 40°-2° 
".^"0.30 

- 10°-°° 

Q-QQ 

- 10°-°° 

Q-OQ 

■ 10°-'° 

Q.QQ 

■ 10°-°° 
"••^"ooo 

■ 10°-°° 
"•'"0.00 

■ 20°-'° 
"•^"0.10 

■ 20°-'° 
"•^"0.10 

■ 20°-"° 
"•^"0.10 

- 30°-'° 

- 30°-2° 
"•-'"0.20 

- 10°-°° 
"•'"0.00 

- 20°-^° 
^-'•^"0.10 

■ 10°-°° 

Q.IQ 

- 50°-^° 

Q.4Q 

■ 40°-*° 

Q-3Q 

- 10°-'° 

Q-OQ 

■ 10°-'° 

Q-OQ 

- 10°-'° 
"•'"0.00 

- 10°-°° 
^-'•'"0.00 

- 20°-2° 

Q-IQ 

- 1 00°-°° 

Q-OQ 

- 10°-'° 
"•'"p.OQ 

- 0.90° 
-0.10^ 

- 10°-'''' 

Q-OQ 

- 10°-'° 

Q-QQ 

■ 10°-°° 

Q.QQ 

- 10°-'° 

Q.OQ 

■ 30°-'° 
"•-"-'0.00 

■ 20°-*° 
"•^"0.10 

- 20°-°° 

Q-IQ 

- 10°-'° 

Q-OQ 

• 10"-'" 
"•'"0.00 

■ 10"-'" 

• Q.IQ 

- 10°°° 



:§8 



'19. 



81. 
'•213. 
70. 
•236. 

196.*„. 
1 0. 



400. 

313.^°- 

625.p; 

200.°- 
123.™- 
400.°-' 
151 

^-'^•52. 

300.°- 

200.°- 

289.2°- 

248.32° 

200.°- 

298 2: 

^^"•139. 

331.™- 
200.°^ 



733 
714 

500.^ 
656.^*- 
719™ 
681.|; 
695.?„,. 
1000.°- 
806.'20 

999. J;" 

1000. °- 

815.!f- 

996 * 
^^"•111. 

701.° 

500.°- 

500.°- 

786.3°; 

877.""- 

412 "2- 

U2. 
998.23, 

990.J«; 



0.90" 



1^ 
0.00 



8:8§ 



430. 



0.95; 



n.no 



C. Gielen et al.: Silicate features in Galactic and extragalactic post-AGB disks., Online Material p 10 



Table A^. Best fit parameters deduced from our full spectral fitting. The abundances of small, medium and large grains of the 
various dust species are given as fractions of the total mass, excluding the dust responsible for the continuum emission. 
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Table A.3. Best fit parameters deduced from our full spectral fitting. The last column gives the continuum flux contribution, listed 
as a percentage of the total integrated flux over the full wavelength range. 
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